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We have obtained EUV spectra between 90 and 255 A of the comets C/2002 
T7 (LINEAR), C/2001 Q4 (NEAT), and C/2004 Q2 (Machholz) near their per- 
ihehon passages in 2004 with the Cosmic Hot Interstellar Plasma Spectrometer 
(CHIPS). We obtained contemporaneous data on Comet NEAT Q4 with the 
Chandra X-ray Observatory ACIS instrument, marking the first simultaneous 
EUV and X-ray spectral observations of a comet. The total CHIPS/EUV ob- 
serving times were 337 ks for Q4, 234 ks for T7, and 483 ks for Machholz and 
for both CHIPS and Chandra we calculate we have captured all the comet flux 
in the instrument field of view. We set upper limits on solar wind charge ex- 
change emission lines of O, C, N, Ne and Fe occurring in the spectral bandpass 
of CHIPS. The spectrum of Q4 obtained with Chandra can be reproduced by 
modeling emission lines of C, N O, Mg, Fe, Si, S, and Ne solar wind ions. The 
measured X-ray emission line intensities are consistent with our predictions from 
a solar wind charge exchange model. The model predictions for the EUV emission 
line intensities arc determined from the intensity ratios of the cascading X-ray 
and EUV photons arising in the charge exchange processes. They are compatible 
with the measured hmits on the intensities of the EUV hnes. For comet Q4, we 
measured a total X-ray fiux of 3.7x10"^^ ergs cm~^ s~^, and derive from model 
predictions a total EUV fiux of 1.5x10"^^ erg cm~^ s~^. The CHIPS observa- 
tions occurred predominantly while the satellite was on the dayside of Earth. 
For much of the observing time, CHIPS performed observations at smaller solar 
angles than it was designed for and EUV emission from the Sun scattered into 
the instrument limited the sensitivity of the EUV measurements. 

Subject headings: comets: general — comets: individual (C/2002 T7 (LINEAR), 
C/2001 Q4 (NEAT), and C/2004 Q2 "Machholz") — ultraviolet: solar system 

1. Introduction 

Soft X-ray emission has been detected from at least 20 comets since 1996. After the 
consideration and rejection of several possible emission mechanisms, a consensus has been 
reached that the primary mechanism is charge exchange collisions between highly charged 
solar wind minor ions and neutral atoms and molecules of the cometary atmospheres, as 
suggested by Cravens (1997). The observations have been reviewed by Cravens (2002), 
Krasnopolsky et al. (2004) and Lisse et al. (2004). The solar wind charge exchange (SWCX) 
process also occurs in the heliospherc where it contributes to the soft X-ray background (Cox 
1988; Cravens 2000; Robertson & Cravens 2003; Pepino et al. 2004; Lallement 2004) and may 
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be responsible for some fraction of the short and long term enhancements seen at 250 eV 
(Snowden et al. 1995, 2004). These recent, interesting observations indicate it is important 
to understand both the microphysics of the charge-exchange process and the phenomenology 
of where, when and how this process occurs within the solar system. 

Theoretical models of the X-ray and EUV emission lines produced in comets by charge 
exchange with the solar wind ions have been constructed by Wcgmann ct al. (1998), Schwadron 
& Cravens (2000), Kharchenko & Dalgarno (2000), Kharchenko & Dalgarno (2001), Kharchenko 
et al. (2003), Beiersdorfer et al. (2003), and Lisse et al. (2004). 

Observations of diffuse SWCX emission other than from comets have been carried out 
using Chandra (Wargclin ct al. 2004) and XMM-Newton (Snowden ct al. 2004). Wargelin 
et al. (2004) observe the optically dark side of the moon and attribute the X-rays they 
detect to charge exchange collisions between solar wind ions and neutral hydrogen in the 
Earth's geocorona. Snowden et al. (2004) ascribe the emission they see as coming from the 
Earth's exosphere or the hehosphere. In both cases, the emission is characterized by strong 
variabihty and correlation with the solar wind particle flux. X-rays generated by SWCX 
have been detected from the atmospheres of Mars and Venus (Murawski & Steinolfson 1996; 
Dennerl 2002; Holmstrom & Kallio 2004; Fok et al. 2004; Gunell et al. 2004). 

In summary, the highly-ionized, heavy ions that stream out in the solar wind are a 
potential source of soft X-ray and EUV emission from within the solar system. Charge 
exchange emission may occur involving neutral atoms and molecules from comets entering the 
inner solar system, planetary exospheres, and neutral atoms in the inter-planetary medium. 
A highly-ionized solar wind ion may give rise to several photons as the ion captures electrons 
sequentially in its flight from the Sun. A detailed ionization and density map of the heavy 
solar wind ions as a function of time and position is required for a complete description 
of the X-ray and EUV emissions in the solar system. Measurements of SWCX emission 
from comets (and other measurements) can provide a useful probe of the density and ion 
composition of the solar wind. 

Many SWCX lines are predicted to appear at EUV energies below the low-energy limit of 
Chandra and XMM-Newton (Kharchenko & Dalgarno 2000, 2001, Kharchenko et al. 2003). 
Some of these EUV hues arise from ions not conclusively identified from emission at higher 
energies. These fines may be detected with CHIPS, a spectrometer sensitive in the EUV with 
resolution of several hundred. The lines appear to be present in the ROSAT PSPC spectrum 
of C/1990 Kl (Levy) ((Dennerl et al. 1997)) which rises steeply at energies below 200 eV. 
Krasnopolsky & Mumma (2001) observed comet C/1996 B2 (Hyakutake) at low resolution 
at EUV wavelengths and claimed detection of lines from O, C, He and Ne. At lower energies, 
Kharchenko & Dalgarno (2001) predict that emission lines from the O VI doublet at 1032 
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and 1038 A should be the brightest FUV hnes SWCX hnes. A sensitive search for EUV and 
FUV emission hnes cometary SWCX was conducted with the Far-Ultraviolet Spectroscopic 
Explorer (FUSE) observatory (Weaver et al. 2002; Feldman 2005). Weaver ct al. (2002) 
report a a tentative detection of 1032 A emission line from comet C/2000 WMl (LINEAR). 
An alternative interpretation of this feature has been suggested by Feldman (2005) from in 
FUSE observations of the brighter comet C/2001 Q4 (NEAT); they attribute this feature 
to H2 fluorescent emission line pumped by radiation from N 111 and O VI ions in the solar 
corona. The non-detection of the brightest FUV O VI hnes from SWCX may be because 
the relatively small field of view of FUSE allowed observation of only about 6000-7000 km 
around the cometary nucleus. For the Q4 FUSE observations, the peak brightness of the 
charge- exchange EUV and X-ray emissions was shifted by more than this amount from the 
comet center, hence the FUSE aperture may have missed it (Feldman 2005). 

The favorable perigee passages of the bright Comets C/2001 Q4 NEAT and C/2002 T7 
LINEAR in May 2004 and Machholz in January 2005 provided an excellent opportunity to 
use the unique low energy spectroscopic capabilities of CHIPS to study SWCX in comets at 
EUV wavelengths. Due to a combination of high solar activity and their close approaches, 
Comets C/2001 Q4 NEAT and C/2002 T7 LINEAR are the brightest X-ray targets available 
since Comet C/Hyakutake in 1996 and Comet 153P/2002 CI (Ikeya-Zhang) in 2002. In this 
paper, we discuss the observational data from CHIPS and Chandra, describe the procedures 
by which we measure the line strengths and discuss the implications for the cometary charge 
exchange process. 



2. Observations and Data Reduction 
2.1. Instrumental 

The CHIPS orbit is approximately circular at 600 km and an inchnation of 94 degrees, 
determined by the primary payload (ICESAT) of the rocket that launched the instrument. 
CHIPS is sensitive to the EUV band of 90 - 260 A, or 47 - 138 eV. With all six slits open, 
the spectrograph field of view is roughly rectangular at 5 by 25 degrees. Any object smaller 
than a few degrees presents a "point source" to CHIPS, and will be visible only through two 
of the slits at a time. The observations presented here were collected in wide-slit mode, with 
a peak resolution for diffuse hght of about 3.4 A FWHM, reaching about 4 A at the edges 
of the bandpass. 

Light from the comet enters the CHIPS spectrograph through two slits, each 1 mm 
by 7 cm. The light diverges onto an 8 cm x 6 cm grazing incidence cylindrical diffraction 
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grating with central groove density of 1800 lines mm~^. The groove density is varied across 
the grating surface to provide nearly stigmatic focusing in the dispersion direction near the 
center of the spectral band. Thus, the spectral resolution near the center of the band is 
limited by the slit width. At the ends of the band, the combined effects of optical aberration 
and defocusing contribute to the net spectral line width. For sources that are significantly 
smaller than the field-of-view of a slit, the resolution is independent of the angular size of 
the comet. 

The grating rulings and baffling reject most non-EUV light and thin-film filters placed 
close to the detector also attenuate out-of-band and scattered light. A diffuse emission 
line spans the full available detector height, typically passing through two different filter 
materials: aluminum, zirconium, or polyimide/boron. Measurements of an emission line 
that has the appropriate ratio of line flux from the two different fllters spanned by the line 
encourages the belief that an actual spectral line has been detected and not some detector 
artifact. 

The CHIPS in-orbit cahbration is discussed in more detail in Hurwitz et al. (2005) 
and additional details of the science instrument and the CHIPSAT satellite are presented in 
Hurwitz et al. (2005), Hurwitz et al. (2003), Janicik et al. (2003), Marckwordt et al. (2003), 
Sirk et al. (2003), and ShoU et al. (2003). 

2.2. Comet Observations and Data Reduction 

A summary of the CHIPS observations of all three comets is presented in Table 1. 
We coadd individual comet observations in different subsets as described below. Individual 
pointings are of ~20 minute duration and owing to the orbital geometry, nearly all obser- 
vations took place while the satellite was on the day side of Earth. The CHIPS mission 
has a nominal constraint of 72 degrees for the minimum borcsight-to-Sun angle permitted 
during observations, both for spacecraft health reasons and to minimize the amount of scat- 
tered solar radiation entering the spectrometer. However, observing NEAT Q4 and LINEAR 
T7 near their time of maximum optical brightness required observations that violated this 
Sun-angle constraint and required the development and use of a new observing mode to 
accommodate the smaller solar angles as well as observe a target moving in the sky. The 
two most important differences for this mode are a slightly reduced pointing accuracy and 
stability since the main Sun sensor did not have the Sun in the fleld of view and a possible 
increase in the amount of scattered solar EUV radiation present in the comet spectra. 

Initial data reduction was carried out using the pipeline that is described in Hurwitz 
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et al. (2005). We apply several filters to the raw data in producing the reduced spectra. 
Data are excluded when the overall detector count-rate exceeds 80 events/sec, indicative of 
high charged-particle backgrounds. Data are also excluded when the detector high voltage 
is reduced (a normal, temporary response to transient high-background periods). A broad 
pulse-height filter is applied by the flight software. Approximately 40% of the telemetered 
events are then excluded by a low-pass pulse-height threshold in ground software, selectively 
reducing background because the low amplitude events are overwhelmingly triggered by 
charged particles rather than by photons. 

Periodically pulsed "stimulation pins" outside the active field of view arc used to register 
the event X,Y coordinates to a common frame, thereby correcting for thermal drifts in 
the plate scale or axis zero-points, and to determine the effective observing time of the 
observation. Distortion is corrected using a preflight pinhole grid map. Small regions of 
known detector "hot spots" are excluded, the event coordinates are rotated so that the 
new X axis closely corresponds to the spectral dispersion direction, and the spectra in each 
filter-half are summed over the active detector height. 

Once the final spectra (one each from the upper and lower detector halves) have been 
summed, these will typically contain counts from several background sources in addition to 
any signal from a comet. The primary backgrounds include particle events that arc generally 
uniform across the detector, geocoronal He 1 584 A photons that are scattered by the gratings 
and penetrate the Al filters, and for these comet observations, some EUV photons from the 
Sun that are multiply-scattered from some structure on the satellite into the spectrometer. 
These backgrounds must be properly accounted for or subtracted before we can assess the 
amount of radiation from the comets that may be present. 

We show in Fig. 1 the summed spectra from all the individual NEAT Q4 observations, 
totaling 329495 s of useful data. The upper curve is the Polyboron/Al side of the detector, 
while the lower trace is the Zr/Al side of the detector, lowered by 200 counts. The wavelength 
ranges of the various filters are indicated, as are the positions of the filter-supporting bars 
for each half. The continuum seen in the Poly/B and Zr filters is primarily from particle 
events, while the Al filters show a significant additional smooth background from scattered 
geocoronal He I 584 A photons. The effective area of the spectrometer as a function of 
wavelength is shown in Sirk et al. (2003). We add together signal from one or both halves of 
the detector at a particular wavelength based on the filter transmission curves to maximize 
signal-to-noise at that wavelength or to exclude artifacts. The ranges included from each 
filter/half are indicated in Fig. 1 . 

To determine the continuum against which to search for emission lines, we use a very 
deep fiight charged-particle spectrum, scaled to match the local counts in the spectrum 
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under analysis. The deep charged-particle spectrum generally contains many more events per 
wavelength interval than does the spectrum under analysis, and thus contributes relatively 
little to the shot noise in the line flux. Differences in the technique by which the scale factor 
is determined - for example, heavy smoothing versus fitting low-order polynomials - result 
in changes of ~10% to the line fiux results. To determine the absolute "fiatness" of the 
detector response, we analyzed a pre- flight photon flat fleld that was histogrammed like the 
flight spectra. Compared to a smooth polynomial, the flat field spectrum showed a pixel-to- 
pixel variance only slightly greater than expected for shot noise, where "pixel" refers to the 
width of the narrowest spectral features. The excess one-sigma variance above shot noise 
was about 0.3% of the total signal. This factor is negligible in all the fimits reported here, 
where shot-noise from charged particles detected is the dominant uncertainty. 

We measure the line fiux of an individual line by adding all the counts that occur 
within one spectral resolution element of the line rest wavelength, and subtracting the counts 
summed from the scaled background spectrum over the same limits. 

The wavelength scale is based on prefiight measurements, offset by a (constant) ~1 A 
determined from the bright geocoronal He II 256.3 A feature, which is present in essentially 
all the observations. Temporal variations in the measured centroid of the He II feature show 
a dispersion of about 0.2 A. Periodic observations of the moon, which provide a weakly 
reflected solar spectrum that includes lines of Fc IX through XII, confirm that the adopted 
wavelength scale provides a good fit near the center of the spectral band, and support the 
view that the relative throughput of the filter panels is as measured prc-fiight. The measured 
fiuxes of both the He II feature and the lunar spectrum are in good agreement with pre-fiight 
expectations, suggesting that the instrument throughput has not declined since laboratory 
cahbration. The pre-fiight laboratory throughput cahbration, tied to an NIST-calibrated 
diode, both at a component level and an end-to-end manner (which agree to within ~10%) 
provide the absolute fiux calibration. 



2.3. Pointing Stability 

The CHIPSat attitude control system (ACS) is described in Janicik et al. (2003). The 
ACS uses inputs from several sensors, include coarse and medium Sun sensors, magnetic sen- 
sors, a lunar sensor, and reaction wheel rotation sensors to calculate the spacecraft attitude. 
During most of the observations reported here, CHIPS was unable to use its most accurate 
source of attitude information, the spacecraft's medium Sun sensor. Instead, the adopted 
observing procedure was to orient the spacecraft such that the Sun was in the Sun sensor 
for the fraction of each orbit when a comet was not being observed, then slew to the comet 
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during the period the comet was visible. When the Sun sensor is used, the typical point- 
ing accuracy of the spacecraft is within ~0.5 degrees of the requested target Janicik et al. 
(2003), and centered to within 1/lOth of the spectrometer field-of-view. For relatively short 
periods even with slews, such as the length of a typical 25 minute comet observation, the 
spacecraft reaction wheels provide pointing orientation and knowledge with negligible degra- 
dation. This is confirmed by the very small attitude corrections that followed reaquisition 
of the Sun by the medium Sun sensor. 

In the CHIPS observations of the comets, we conclude that essentially all of the EUV- 
emitting region of the comet is contained within the CHIPS entrance apertures during the 
observations. The field of view of the two central slits of CHIPS is approximately five 
degrees square. Comet positions were generated every 10 minutes during the period they 
were observed by CHIPS using the Minor Plant and Comet Ephemeris Service^, where the 
requested pointing corresponded to the mid-point in time of each observation. Parallax and 
the motion of the comet during an orbit's pointing are negligible compared with the CHIPS 
field of view. Our re-mapped image of Q4 on the Chandra ACIS-S detector shows the X-ray 
emitting region to be no larger than 10 x 10 arcmin. Given the attitude information and 
the size of the X-ray emitting region, we believe that the comets were positioned in the 
slit entrance aperture when the spacecraft reported it was pointed at the requested comet 
position and that the emitting region was fully contained within the CHIPS field of view. 
Hence, for the observations reported here, we assume no dilution factor for comet EUV 
emission. The roll angle for observations was chosen such that the CHIPS entrance slits 
were perpendicular to the great circle that contained both the comet and the Sun. 



2.4. Background/Noise Sources 

The majority of counts registered in the CHIPS detectors are from charged particles 
entering the instrument. Although some discrimination against particle events based on 
pulse height takes place as described above, the time-varying nature of this background 
makes it difficult to determine the presence of any continuum emission arising from photons. 

The net effect of the backgrounds and data reduction techniques is that CHIPS is primarily 
sensitive to spectral lines, with a limiting sensitivity of a few line units or LU [photons cm~^ 
g-i s];-i] fQ]; lY^Q observations described here. 

CHIPS consistently detects at high significance He II emission arising in the Earth's 
geocorona at 256.3 A and to a lesser degree at 243.0 A. In addition, we observe radiation 



^http://cfa- www.harvard.edu/iau/MPEph/MPEph.html 
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that appears to originate in the solar corona and is multiply scattered from surfaces on the 
spacecraft into the instrument. Although this background is very faint compared with the 
interstellar medium emission lines from Fe IX - XII that CHIPS was designed to observe, we 
have studied it in depth because it is a signficant background to these measurements. The 
EUV line emission from the interstellar medium reported by Hurwitz et al. (2005) based 
on one year of CHIPS observations are very faint and do not pose a significant background 
for the comet observations. A detailed study of the instrument-scattered background shows 
that the lines are never seen in co-added nighttime-only observations; the spectrum appears 
very similar to the full moon spectrum; the strength of the lines does not correlate with 
zenith angle; the lines show no ram-angle dependence; the lines appear reproducibly in only 
a relatively small number of spacecraft orientations relative to the Sun; the lines appear more 
often, but not always, when CHIPS points closer to the Sun than during normal observations, 
but there is not a significant anti-correlation with solar-target angle; the line strengths show 
no correlation with the Ultra Low Energy Isotope Spectrometer (ULEIS) (Mason ct al. 1998) 
measurements of Fe IX, X and XI in the solar wind during the time of the comet observations; 
and the line strengths during comet observations do not correlate with measurements of solar 
disk irradiance in the 170 - 230 A by the Solar EUV Experiment on TIMED (SEE) (Woods et 
al. 1998). We conclude that the source of these background lines is faintly-scattered radiation 
from the solar corona. Owing to its appearance preferentially at smaller solar angles, this 
background is therefore of more importance for the CHIPS comet observations than during 
normal sky survey observations, which are made at larger solar angles. 

Fig. 2 shows the co-added lunar spectrum observed by CHIPS in narrow slit mode from 
several months in 2003, and the co-added spectrum from a number of observations when 
the background lines were strongest. The similarity of the spectra is evident, although wc 
would not expect identical spectra because the solar spectrum itself varies and the scattered 
spectrum was taken in widc-sIit mode, which is more sensitive to the scattered radiation 
and has lower spectral resolution. We can nonetheless use either the moon spectrum or 
the co-added bright background times as a template to indicate the possible presence of 
solar-scattered hues in the following way. We developed this indicator by first subtracting 
the continuum level from the scattered light spectrum after fitting it with a polynomial. 
We then similarly fit and subtract the continuum from another CHIPS observation. When 
the two continuum-subtracted spectra arc multiplied together in the region of the lines, the 
template acts as a matched filter sensitive to the solar lines. We call the exposure time- 
normalized product of the two spectra integrated over the spectral region 170 - 230 A the 
"line flux indicator," or LFI. We find LFI to be a sensitive indicator of the presence of 
scattered solar hght in co-added observations and can use it to reject data that contains this 
background where required. The third trace in Fig. 2 is co-added daytime data selected for 
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low LFI and shows essentially no solar contamination. However, in most cases we measure 
only upper limits to the spectral lines expected from charge exchange reactions, so that 
the presence of a possible solar background is of less concern than it would be were CX 
lines detected. Additionally, the CX lines we are trying to detect form an only partially 
overlapping set with the lines emerging from the solar corona. 

We show in Fig. 3 the relative LFI value for each day of the Q4 and T7 CHIPS ob- 
servations. We also show the scaled magnitude of the solar 170 - 230 A flux of the Sun as 
measured by SEE, which does not correlate with the LFI values. Finally, we also show a 
model visual brightness of the two comets, where the variation is based only on the inverse- 
squared distances from the Earth to the comet and the comet to the Sun. The latter curves 
indicate that we observed Q4 through its peak brightness, while the T7 observations were 
only possible as the brightness dechned, since the comet was too close to the Sun in the sky 
to observe prior to that time. Owing to the sensitivity level of the CHIPS data, it is not 
useful to subdivide the data into different time slices to investigate possible time variation 
of the emission. 



2.5. Chandra observations of C/2001 Q4 (NEAT) 

C/2001 Q4 (NEAT) was observed by Chandra May 12, 2004 using director's discre- 
tionary time. The total approved observation of 10 ks was divided into 3 pointings of 
approximately equal length. During each observation the telescope was held fixed and the 
comet allowed to drift through the field of view. The comet was centered on the S3 chip 
which has the best low energy response. The three resultant standard level 2 event lists 
were simultaneously concatenated and reprojected into a fixed "comet-centered" frame of 
reference using the CIAO^ tool " sso -freeze" . The final reconstructed image provides good 
spatial ('^0.5" /pixel) and spectral information (a ~50 eV). 

A source spectrum was extracted from a circular aperture with radius of 700 pixels (5.8') 
centered on the peak of the comet emission and the background spectra was extracted from 
the area outside the source aperture. We examined how varying the aperture radius changed 
the comet and background spectrum. The aperture radius we chose included essentially all 
of the emission from the comet; outside of this radius we expect that hehospheric SWCX 
will be a uniform background on this scale. The background spectrum normalized by area 
was subtracted from the in-aperture spectrum to produce the Q4 comet spectrum. Work by 
Lisse et al. (1996) and Wegmann et al. (2004) shows the CX X-ray emission from a comet 
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as active as Q4 may be offset from the nucleus by 30,000 - 40,000 km. At the distance 
of 0.6 au where Q4 was observed with Chandra, this separation translates to an angle of 
~1 arcminute. Owing to the diffuse nature of the emission and the positional uncertainty 
induced by the photon remapping it is not clear whether we can discern a potential offset 
of this magnitude from the center of the nucleus, but we have verified that the spectrum 
resulting from our chosen aperture radius is insensitive to a potential offset of this size. 
A detailed discussion of the Q4 X-ray emission morphology will be covered in Lisse et al. 
(2007). ACTS response matrices to model the instrument effective area and energy dependent 
sensitivity were created with the standard CIAO tools and the resulting spectra were fit using 
XSPEC. 



3. Results and Discussion 

Fig. 1 shows the comet C/2001 Q4 (NEAT) spectrum from 50 - 270 A measured with 
CHIPS. For each potentially strong charge exchange spectral fine indicated in the figure, we 
measure the brightness using the spectra from one or both halves of the CHIPS detector. 
We have run our line fiux measuring algorithm for each comet and the results of these 
measurements are shown in Table 2. Wc include two measurements using Q4 data: the 
"no LFI" set refers to the period May 9 - 19, 2004 corresponding to a time when the LFI 
measurement indicates minimal contamination from scattered solar EUV radiation. We also 
present the limits from a spectrum consisting of the co-added spectra from all three comets, 
which contains the lowest relative shot noise background. We measure positive fiux at the 
wavelengths of a number of potential charge exchange lines, albeit at low significance. 

Fe lines at wavelengths 170 - 220 A are detected in several of the spectra at significances 
ranging to 5 a. Since it is likely that much or all of this emission is solar rather than 
cometary, we restrict our consideration to possible charge exchange lines shortward of 170 
A, although line measurements that may have a solar component provide upper limits to 
the cometary emissions. We estimate that the lines for which we report limits in the table 
contribute about 90% of the fiux from charge exchange if we ignore the iron lines. 

The least contaminated spectrum is that obtained from the NEAT Q4 data selected for 
low LFI and the most significant line is the O VI line at 129.91 A, at almost 2a, arising from 
the ls4d - ls2p ^P transition. 

In Fig. 4 we show the results, binned to the CHIPS instrumental resolution, of the line 
fiux measuring algorithm run on a 1 A wavelength grid on the combined data set. When 
using this technique a positive actual line fiux suppresses the neighboring continuum. Hence 
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positive flux should be read from the graph axis, not the adjacent continuum. Tentative 

hne identifications of features associated with positive flux measurements are indicated. The 
feature at 91.1 A coincides with the transition ls4p '^P - ls2s '^S of O VII, and the feature 
at 128.5 A to the ls3d - ls2p transition. O VI has hues at 116.4 A, 129.9 Al32.5 A 
and 150.1 A, corresponding respectively to the 4p-2s, 4d-2p, 4s-2p and 3p-2s transitions. 

A search for line detections can also be conducted using theoretical predictions for the 
wavelength and flux ratio of several lines from the same ion. We employed the nine brightest 
lines in the CHIPS band of the O VI ion given by the model described in the following section. 
The line profllcs were taken to be Gaussian with widths equal to the wavelength-dependent 
resolution of CHIPS. The free parameters in the flt were a quintic polynomial to find the 
background level and an overall amplitude for the O VI lines, with the individual line ratios 
fixed. The results when this technique was apphed to the NEAT Q4 "no LFI" spectrum are 
that the line amplitudes are consistent with zero. We carried out the same procedure for the 
co-added comet T7 spectra and again the results were consistent with zero O VI emission. 

4. Modeling the Solctr Wind Charge Exchange Emission 

The NEAT Q4 count spectrum obtained from CHANDRA is shown in Fig. 5. It 
contains emission fines seen in other comets (Lisse 2001; Krasnopolsky et al. 2002). We first 
use XSPEC to fit the spectrum with eight emission fines whose wavelengths were allowed 
to vary. The best fit was achieved with 8 lines at energies in cV of 245, 320, 403, 496, 570, 
655, 830 and 952, with corresponding photon fluxes in units of 10"*^ photons cm~^ s~^ of 
respectively 3.68, 26.5, 6.57, 2.59, 5.93, 1.05, 0.29, and 0.11. The resulting per degree 
of freedom was 0.91. The strong lines at 320, 403 and 570 eV correspond to transitions of 
C V, N VI and O VII. The lines near 245 eV may be from Si IX, Mg IX or Mg X, and 
a feature at 952 eV found also in the spectra of Comets McNaught-Hartley and Linear 4 
(Krasnopolsky et al. 2002, 2004) may be from Ne IX. There is an indication of a line near 
835 eV but we hesitate to attribute it to the 5p-ls resonance fine of O VIII because we do 
not find the corresponding 2p-ls fine at 653 eV which we predict should be a factor of six 
more intense. Because of the low resolution of the instrument, fits with a similar ^ figure of 
merit are also be obtained by a combination of emission lines and thermal bremsstrahlung 
at a temperature of about 0.23 keV, as in Kharchenko & Dalgarno (2000). 

We have extended a model developed by Kharchenko &: Dalgarno (2000), Kharchenko 
& Dalgarno (2001), Kharchenko et al. (2003) and Pepino et al. (2004) to compute emission 
spectra over the full EUV and X-ray spectral range. In the model we use electron capture 
cross sections for positive ions colliding with neutral H2O molecules that reproduce the ex- 
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perimental laboratory X-ray spectra when available and theoretical estimates when they are 
not (Rigazio et al. 2002). When applied to the measured spectra the method takes into ac- 
count empirically-determined processes such as multiple capture followed by autoionization. 
The capture cross sections determine the entry rates into the individual excited states of the 
ion created in the charge exchange. The excited states decay in a radiative cascade which 
produces photons at X-ray, EUV and UV wavelengths. The relative intensities of the emis- 
sion lines occurring in the cascade from a given initial level depend only on the branching 
ratios of the transition probabilities. Thus the EUV spectrum is determined by the X-ray 
spectrum and the comet spectra depend on and provide a measure of the solar wind ion 
composition. 



4.1. X-ray emission lines 

We calculated the X-ray emission line intensities for the ion composition listed in Table 
3. It is based on the ion composition of Schwadron & Cravens (2000) for a fast solar wind, 
supplemented by the presence of a fraction 0.001% of highly charged Fe ions, Fe XV-Fe XX, 
that appear to be necessary to reproduce the observed spectrum between 0.8 and 0.9 keV. 
An enhanced flux of highly-charged iron ions may be a consequence of fast coronal mass 
ejection events (Lepri & Zurbuchen 2004; Zurbuchen et al. 2004). 

A comparison of the observations in counts per second of Comet NEAT Q4 and the 
XSPEC and SWCX models at photon energies between 0.3 kcV and 1.0 kcV is presented in 
Fig. 5. The CHANDRA data with the error bars indicating 1 sigma uncertainties are shown 
as solid circles. The dashed line is the XSPEC model fit with 8 emission lines folded through 
the instrument response. The solid line is the SWCX prediction at a spectral resolution of 
75 eV. In Table 4 we list the energies and wavelengths and the predicted relative intensities 
of the principal X-ray and EUV emission hues of the SWCX model. The agreement of the 
measurements of the X-ray lines with the models is good, but because of the low spectral 
resolution, we cannot exclude a significant contribution from bremsstrahlung from XSPEC 
modehng of the X-ray lines alone. 



4.2. EUV emission lines 

The photon fluxes measured for comet NEAT Q4 with CHIPS at energies below 150 eV 
and with CHANDRA at energies between 250 eV and 1.0 keV are compared in Fig.6 with 
the predictions of the SWCX model. A scahng factor has been applied to bring the predicted 
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total flux of the X-ray emission with energies greater than 250 eV into agreement with the 
measured value of 3.7x10"^^ ergs cm~^ s~^. Table 4 can then be used to yield the absolute 
intensities shown in Fig. 6. The corresponding total flux of the EUV lines between 47 eV and 
140 eV is 1.5 xlO^-*^^ ergs cm~^ s"-*^, distributed amongst 181 lines of C,N,0, Ne and Fe. The 
theoretical EUV intensities fall well below the measurements except near 95 cV, a region in 
which several lines of OVI and OVII appear. The most interesting CHIPS measurement of 
the SWCX hnes is the 2 sigma detection in the NEAT Q4 "no-LFI" spectrum of the line of 
OVI at 95.5 eV or 130 A. Discrepancies at lower energies are affected by scattered radiation 
from the Sun. We conclude that the SWCX model EUV predictions are consistent with 
the CHIPS data, but the instrument lacks the sensitivity needed to systematically test the 
theoretical predictions. In that CHIPS was not designed for measurements of comets, or point 
sources, the limited agreement supports the SWCX model and points out some additional 
utility of the CHIPS measurements. The model also predicts that the 10 brightest EUV lines 
in the CHIPS band contribute 65% of all the EUV photons, and 88% of the EUV photon 
flux is contained in the brightest 26 lines. 

The accumulated yield is useful in providing an estimate of the photon energy distribu- 
tion from comets. It is defined as the total number of photons produced by charge exchange 
in collisions of the solar wind with comets with energies lying between a minimum energy 
,equal to 47 eV for CHIPS, and a maximum energy E. The accumulated yield is dimension- 
less. Multiplied by the product of the solar wind flux and ion density it yields the volume 
emissivity. Fig. 7 shows the calculated accumulated yield for comet NEAT Q4 as a function 
of photon energy E for a minimum energy of 47 eV. Multiplying it by a factor of 26 yields 
the accumulated photon flux in units of cm~^ s~^. The small values reflect the low fractional 
abundances of the heavy solar wind ions. 



5. Conclusions 

We have measured EUV emission from three comets and provided the flrst-ever simul- 
taneous EUV and X-ray measurement and modehng of cometary SWCX emission. The 
X-ray data for comet Q4 show clear emission lines that are satisfactorily interpreted as solar 
wind charge exchange lines from ions expected in the solar wind. We have used the same 
model to identify and predict the strengths of EUV lines. The measured EUV line strengths 
are consistent with the theoretical model that reproduces the X-ray line intensities. The 
CHIPS data and model results show that the EUV lines from all three comets are quite 
faint. Nonetheless, the high spectral resolution measurements of the EUV lines have now set 
useful limits on many emission lines in this band that may serve as a guide for future EUV 
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observations of comets. 
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Table 1. Summary of CHIPS comet observations 



Comet 


Date Range Obs. 


Exp. Time 


Exp. Time 


Sun Angle 


Sun-Comet 


Earth-Comet 


V mag. Range 


0.25-1 KeV Flux 






Day (s) 


Night (s) 


Range (deg) 


Distance Range 


Distance Range 




10~^^ ergs cm~^ s~^ 




(CHIPS) 








(AU) 


(AU) 






C/2001 Q4 (NEAT) 


2004^04^21 07:35 - 2004^05-21 07:39 


331096 


5832 


78.1 - 69.0 


1.05 - 0.97 


0.58 - 0.56 


2.6 - 2.1 


3.7 


C/2002 T7 (LINEAR) 


2004^05-02 11:45 - 2004^05-31 03:41 


219477 


14415 


40.0 - 72.8 


0.65 - 1.0 


0.76 - 0.55 


1.5 - 2.7 




C/2004 Q2 (Machholz) 


2005-01-06 23:44 - 2005-01-28 22:58 


441466 


41158 


130.3 - 106.7 


1.24 - 1.21 


0.35 - 0.47 


4.1 - 4.7 





Column 9 from Chandra observations of May 9, 2004 



I 

I— ' 

00 

I 



Table 2. Charge Exchange Line Flux Limits 



C/2001 Q4 no LFI C/2001 Q4 C/2002 T7 C/2004 Q2 All Combined 



Wavelength (A) Ion 


Limit 




(T 


Sign. 


Limit 




(7 


Sign. 


Limit 




(T 


Sign. 


Limit 




(T 


Sign. 


Limit 




(T 


Sign. 


98.37 


Ne VIII 


-0.22 


1, 


.63 


0, 


.00 


-0.87 


1, 


.03 


0, 


,00 


-2.58 


1, 


.23 


0, 


.00 


-1.00 


0, 


,89 


0, 


.00 


-1.30 


0, 


.63 


0, 


,00 


115.86 


O VI 


-0.51 


0, 


.77 





.00 


0.26 





.48 


0, 


,53 


1.01 





.58 


1, 


.74 


0.19 


0, 


,42 


0, 


,45 


0.39 





.30 


1, 


,31 


120.33 


O VII 


0.20 


0, 


,76 


0, 


,27 


0.14 


0, 


,47 


0, 


,29 


-0.83 


0, 


,56 


0, 


,00 


-0.62 


0, 


,41 


0, 


,00 


-0.42 





,29 


0, 


,00 


128.48 


O VII 


0.26 


1, 


,07 


0, 


,25 


0.06 


0, 


,66 


0, 


,09 


0.10 


0, 


,78 


0, 


,12 


0.64 


0, 


,56 


1, 


,13 


0.33 





,40 


0, 


,84 


129.91t 


O VI 


2.22 


1, 


,16 


1 


,92 


0.52 


0, 


,71 


0, 


,73 


-0.34 





,83 


0, 


,00 


0.56 


0, 


,60 


0, 


,93 


0.35 





,43 


0, 


,82 


132.35t 


O VI 


0.56 


1, 


,24 


0, 


.45 


0.83 


0, 


.76 


1, 


,10 


0.14 


0, 


.89 


0, 


,15 


-0.29 


0, 


,64 


0, 


,00 


0.17 


0, 


.46 


0, 


,36 


133.92 


N VII 


-0.04 


1, 


,29 





.00 


1.07 


0, 


.80 


1, 


,34 


0.40 


0, 


.94 


0, 


.42 


-0.69 


0, 


,67 


0, 


.00 


0.11 


0, 


.48 


0, 


,24 


135.02 


C VI 


-0.27 


1. 


,31 


0, 


.00 


0.27 


0, 


.81 


0, 


,34 


0.34 


0, 


.95 


0, 


.36 


-0.08 


0, 


,68 


0, 


.00 


0.12 


0, 


.48 


0, 


,26 


150.21 


O VI 


-1.84 


1, 


,60 





.00 


-0.21 





.97 


0, 


,00 


-1.09 


1, 


.14 


0, 


,00 


0.15 


0, 


,82 


0, 


,19 


-0.24 





.58 


0, 


,00 


173.29 


O VI 


1.84 


2, 


,02 


0, 


,91 


5.16 


1, 


,26 


4, 


,10 


3.48 


1, 


,57 


2, 


,22 


1.85 


0, 


,87 


2, 


,13 


3.27 


0, 


,73 


4, 


,49 


182.28 


C VI 


0.64 


2, 


,56 


0, 


.25 


2.76 


1, 


.63 


1, 


,69 


6.13 


2, 


.11 


2, 


,91 


-0.69 


1, 


,02 


0, 


,00 


1.90 


0, 


.95 


2, 


,00 


184.19 


O VI 


-2.95 


2. 


,56 


0, 


.00 


-1.96 


1, 


.63 


0, 


,00 


-7.77 


2, 


.10 


0, 


.00 


0.55 


1, 


,03 


0, 


.54 


-2.06 


0, 


.95 


0, 


,00 


186.74 


C V 


2.15 


2, 


,60 





.83 


-0.47 


1, 


.65 


0, 


,00 


-0.47 


2, 


.14 


0, 


,00 


1.59 


1, 


,04 


1, 


.53 


0.48 





.97 


0, 


,50 


197.05 


C V 


0.11 


2, 


,80 





.04 


6.09 


1, 


.80 


3, 


,37 


11.65 


2, 


.36 


4, 


,93 


0.50 


1, 


,12 


0, 


,45 


4.72 


1, 


.05 


4, 


,48 


209.44 


N V 


2.15 


3, 


,20 


0, 


,67 


-2.88 


2, 


,03 


0, 


,00 


-3.61 


2, 


,61 


0, 


,00 


1.81 


1, 


,27 


1, 


,43 


-0.87 


1 


,18 


0, 


,00 


227.20 


C V 


4.84 


3, 


,14 


1, 


.54 


3.42 


1, 


.96 


1, 


,75 


5.15 


2, 


.43 


2, 


,12 


-1.01 


1, 


,23 


0, 


,00 


1.75 


1, 


.13 


1, 


,55 


248.75 


C V 


8.32 


4, 


,38 


1, 


.90 


2.06 


2, 


.74 


0, 


,75 


7.36 


3, 


.40 


2, 


.16 


-6.14 


1, 


,76 


0, 


.00 


-0.57 


1, 


.60 


0, 


,00 



Table 2 notes: Limit and local c values in LU. Sign. = Limit/cr. To find a limit in ph cm ^ s ^ multiply the value in LU by 0.0106. 

flndicated lines and those with A > 170 A axe subject to contamination from solar Fe VIII - Fe XII emission 



Table 3. Solar Wind Ion Composition 



Atom I Ionization State 





5+ 


6+ 


7+ 


8+ 


9+ 


10+ 


11+ 


12+ 


13+ - 20+ 


c 


0.058 


0.004 
















N 


0.0074 


0.0012 


0.0009 



















0.053 


0.004 














Ne 








0.006 


0.00006 










Mg 










0.0025 


0.0017 








Si 








0.0015 


0.0014 


0.00007 








s 










0.0004 










Fe 










0.0015 


0.015 


0.0001 


0.0003 


0.001 



Table 3 notes: Fractional abundance by number density for each ion used for X-ray and EUV model spectrum. Blank 

columns indicate zero abundance was used for the model spectrum. 
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Table 4. Relative Intensities of Emission Lines 



Energy (eV) Rel. Intensity Ion Transition 



47.65 


0.19 


C V 


ls3s - 


ls2p =^P 


49.85 


0.34 


C V 


ls3d - 


■ ls2p 


54.58 


0.30 


C V 


ls3p - 


ls2s 


71.56 


0.50 


VI 


3d - 2p 




82.56 


0.20 


VI 


3p - 2s 




93.70 


0.31 


VI 


4s - 2p 




95.43 


0.19 


VI 


4d - 2p 




107.03 


0.41 


VI 


4p - 2s 




298.96 


1^ 


C V 


ls2s - 


Is^ 


304.25 


0.07 


C V 


ls2p - 


Is^ 


307.91 


0.16 


C V 


ls2p ip - 


Is^ 


354.52 


0.07 


C V 


ls3p ip - 


1S2 IS 


367.35 


0.11 


C VI 


2p - Is 




370.92 


0.03 


C V 


ls4p ip - 


ls2 IS 


419.80 


0.03 


N VI 


ls2s - 


1S2 ^S 


435.38 


0.03 


C VI 


3p - Is 




500.01 


0.03 


N VII 


2p - Is 




560.86 


0.14 


VII 


ls2s - 


Is^ ^S 



Table 4 lists the transitions with the highest emission by photon number in the model 
spectrum. The electronic configurations shown in the table for Li-like and H-like ions 
represent transitions for all states of electronic doublets. 

^The intensity of the C V line corresponds to a photon 
flux of 3.3x10"^ photons cm^'^s^^. 
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Fig. 1. — Complete CHIPS spectrum after pipeline processing using all data from C/2001 
Q4 (NEAT) on the two detector halves, shown in counts. The lower spectrum has been 
displaced by 200 counts for clarity. The filter bandpasses and edges, location of known solar 
and geocoronal hues, and detector/hght artifacts that show up in all CHIPS spectra are 
marked. Vertical marks show the location of charge exchange lines predicted in this band 
from comets. Our flux detection limits for a particular line use one or both halves (indicated 
by lines at the top of the figure, P/Al top, Zr/Al bottom) to optimize the signal-to-noise at 
that wavelength. 
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Fig. 2. — CHIPS daysidc comet and lunar spectra through Al (polyimide/Boron detector 
half) filter and standard pipeline processing. The spectra are not divided by the filter 
response. The upper trace shows a set of co-added CHIPS observations selected when the 
dayside background lines are brightest. The middle trace shows a set of dayside spectra 
chosen with low LFI, showing no solar contamination. The lowest trace shows the co-added 
full moon spectrum taken for calibration purposes by CHIPS in narrow slit mode. The lines 
visible in the lunar spectrum arise from the solar corona being reflected by the moon. The 
similarity in the spectral lines points to an origin in the solar corona for the dayside lines in 
the non-lunar spectra as well. 
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Fig. 3. — The figure shows several quantities related to the faint, variable background versus 
day of year, 2004 in order to investigate correlations with the background. The crosses 
indicate relative LFI (see text) strength by day for C/2001 Q4 (NEAT) observations and 
the asterisks for C/2002 T7 (LINEAR) observations. The last six Q4 points are designated 
"no LFI" points. The solid line is a scaled and offset measure of the SEE 170 - 230 A 
flux. The relative observing geometry appears to be a larger determinant of LFI magnitude 
than the solar EUV flux as measured by SEE. We also show, to indicate the relative comet 
brightness where our observations were made, a geometry-only prediction (based on Horizons 
ephemeris) of comet flux based on distance from Earth and the Sun. The dot-dash line is 
for Q4 and the dashed line is for T7. 
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Fig. 4. — All three comet spectra summed, then continuum subtracted. This procedure can 
improve the visibility of strong emission lines, but is subject to uncertainties in the continuum 
level placement. Spectra and errors shown at 4 A resolution, close to the actual instrument 
resolution and the error bars are 1 sigma. Some tentative spectral line identifications are 
indicated. Longward of 170 A, solar Fe fines dominate the spectrum. 
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Fig. 5. — Two models of X-ray emission lines are compared with the Chandra spectrum 
for C/2001 Q4 (NEAT). The solid points are the number counts from Chandra data. The 
dashed curve is the Chandra data interpreted using an XSPEC fit with 8 lines between 245 
eV and 952 eV folded through the instrument response. The solid line shows the SWCX 
model predictions at a resolution of 72 eV. While both models yield and acceptable fit from a 
standpoint, only the SWCX model is based on substantially all of the underlying physics 
of the emission process. 




Fig. 6. — The intensity of the EUV and X-ray emission induced by the solar wind ions 
interacting with the atmosphere of Comet C/2001 Q4 (NEAT). The sohd curve shows the 
prediction of the CX model for the photon flux in the energy interval of 0.047 - 1.1 keV. 
Absolute intensities of the theoretical CX spectra are normalized to the total number of 
X-ray photons observed with the Chandra telescope in the energy interval 0.25-1 keV. The 
triangles with error bars represent Chandra observations and the squares show the result of 
observations with CHIPS. Theoretical X-ray spectra are computed for a FWHM of 70 eV, 
which corresponds to the Chandra observations. The resolution of the EUV spectral lines 
between 47 and 138 eV has been taken to be equal to the CHIPS resolution of 4 A. The dot- 
dash line is the model fit to the unsubtracted comet X-ray spectrum including both comet 
and heliospheric charge exchange components; it shows a predicted absolute upper limit to 
the SWCX emission. The EUV part of this curve is displayed with the same resolution as 
the X-ray portion. 
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Fig. 7. — The predicted Accumulated Yield of EUV and X-ray photons for a single charge 
exchange collision of the heavy solar wind ions calculated from the SWCX model. The 
figure shows the average number of photons induced with energies between 47 eV and E[eV] 
in a single CX collision of the solar wind ions with molecules of a cometary gas. The 
averaging procedure takes into account the relative abundances of the solar wind protons, 
alpha particles and heavy ions, and individual CX cross sections. 



